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Topological phases of matter are described universally by topological field theories in the same 
way that symmetry-breaking phases of matter are described by Landau-Ginzburg field theories. 
We propose that topological insulators in two and three dimensions are described by a version 
of abelian BF theory. For the two-dimensional topological insulator or quantum spin Hall state, 
this description is essentially equivalent to a pair of Chern-Simons theories, consistent with the 
realization of this phase as paired integer quantum Hall effect states. The BF description can be 
, motivated from the local excitations produced when a tt fiux is threaded through this state. For the 

OA ' three-dimensional topological insulator, the BF description is less obvious but quite versatile: it 

^ ' contains a gapless surface Dirac fermion when time-reversal-symmetry is preserved and yields "axion 

^ , electrodynamics", i.e., an electromagnetic E ■ B term, when time-reversal symmetry is broken and 

■ the surfaces are gapped. Just as changing the coefficients and charges of 2D Chern-Simons theory 

' allows one to obtain fractional quantum Hall states starting from integer states, BF theory could 

' also describe (at a macroscopic level) fractional 3D topological insulators with fractional statistics 

of point-like and line-like objects. 
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I. INTRODUCTION 
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Most ordered phases of condensed matter can be classified using the concept of symmetry breaking, which leads to a 
I ■ theoretical description in terms of Landau-Ginzburg field theories. Such theories have had remarkable successes, such 
' ^ ] as the quantitative explanation of critical phenomena. Topological phases such as the quantum Hall effect require a 
different type of description in terms of "topological field theories" . These theories have the remarkable property that 
they are independent of the spacetime metric. The most famous example is the description of quantum Hall phases 
by Chern-Simons theory, which we review in a moment. As shown by Wen, that theory leads to detailed predictions 
about edge tunneling experiments on fractional quantum Hall samples, which to some degree have been verified at 
least qualitatively. 

Recently new topological phases of electrons with time-reversal symmetry have been discovered. These "topological 
insulators"-'- result from spin-orbit coupling and exist both in two'^"^ and three^^^ dimensions, unlike the quantum 
Hall effect, which is essentially two-dimensional. These phases are insulating in bulk but support conducting edge 
' or surface states. Hence one picture of a topological insulator is as a bulk insulator that necessarily has conducting 
edge or surface states; more recently there has been considerable interest in the electromagnetic responsei^iii when a 
weak time-reversal-breaking perturbation gaps the surface state of the three-dimensional (3D) topological insulator. 
Our goal in this work is to find an effective field theory description of topological insulator phases in the same spirit 
as the Chern-Simons effective theory of the quantum Hall effect. 

In order to understand why we do not consider the electromagnetic response in the 3D topological insulator to 
be a topological field theory in the sense of Chern-Simons theory, it is useful to review the quantum Hall case. In 
the quantum Hall effect, the Chern-Simons theory appears in two forms. We assume for now a single-component 
^ . theory, such as describes the Laughlin states at v — 1/k. In terms of an internal gauge field and an integer fc, the 
Chern-Simons Lagrangian density is 

Lcs - ^e'^'^a^.a.flA. (1) 
Upon coupling this term to electromagnetic fields by adding the term 



-J^A^, = ^e^-^d,ax (2) 



27r 

integrating out the internal gauge field generates a term for the electromagnetic gauge field that has the same 
Chern-Simons form as the internal gauge theory, but with coefficient rather than (Note that A^ does not 

include the background magnetic field generating the quantum Hall state.) 

So for the quantum Hall effect, both the internal theory and the electromagnetic term have the same Chern-Simons 
form. The Chern-Simons theory on a closed two-dimensional manifold has a set of zero-energy states, where g 
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is the genus of the manifold, and an infinite gap to other states. The quantum Hall response is generated by the 
ground state's response to electromagnetic perturbations in similar fashion to the superflow in a superconductor. 
On a manifold with boundary, the bulk Chern-Simons theory alone is not gauge-invariant, because under a gauge 
transformation, ((IJ changes by a total derivative. This total derivative leads to the topological part of the gapless 
edge theory; combining it with non-universal interactions from the edge gives the chiral Luttinger liquid description 
of edge states. 

Our general goal is to carry through a similar program for topological insulators. We note that many features 
of three-dimensional topological insulators with gapped surfaces have been understood in terms of a total-derivative 
term for electromagnetic fields known as "axion electrodynamics" ^'^i^^ : 

This is analogous to the electromagnetic Chern-Simons term for in the quantum Hall effect. It is "topological" 
in two senses: it is a total derivative of the electromagnetic Chern-Simons term, indicating that the surfaces where 9 
changes support a half-integer quantum Hall effect. Its microscopic origin in a crystal is definitely topological as it 
involves the Chern-Simons form of the Berry connection of the Bloch electronsiSiiiii^iii. 

We will argue that the natural description of the internal gauge theory takes a different form, which is one important 
difference between the quantum Hall and topological insulator cases. Note that one difference between the axion 
coupling and the Chern-Simons term is that the axion term is gauge-invariant, while the Chern-Simons term 
is only gauge-invariant up to a boundary term. So any theory based purely on the axion coupling for an internal 
gauge field cannot generate a dynamical boundary from the requirement of gauge invariance in the way that the 
Chern-Simons theory does. 

Our approach leads to a version of BF topological field theory in two and three dimensions as the effective description 
of topological insulators. This theory, which we introduce starting with the two-dimensional (2D) case in the following 
section, preserves time-reversal and parity symmetries, unlike the Chern-Simons theory; in particular It also can exist 
in three dimensions (3D), again unlike the Chern-Simons theory. In two dimensions, it is well known that the relevant 
BF theory can be recast as two copies of Chern-Simons theory^^'^^, but we will argue that the BF picture captures 
better the important properties of the 2D topological insulator, in particular its intrinsic time-reversal invariance. 

We obtain the 2D BF theory heuristically in the following section starting from the distinctive response of a 2D 
topological insulator to an electromagnetic tt flux (i.e., insertion of a solenoid with flux hc/2e)r^—. In 3D, the 
BF theory involves one vector gauge field and one rank- two tensor field These fields are associated with the 
distinctive response of a 3D topological insulator to a fiux linc^^ , which is closely related to the response of the "weak" 
topological insulator to a dislocation2S.. So in both 2D and 3D we can understand BF theory from the attachment of 
nontrivial responses to electromagnetic flux. 

Section HI introduces some general features of BF theories of relevance to topological insulators, and Section 
IV provides more details for the 2D case. Section V contains most of the main results of this paper on the 3D 
case. One type of electromagnetic coupling is present even with time-reversal symmetry, while breaking time-reversal 
symmetry at the boundary is shown to lead to the EB term after integrating out the internal fields of BF theory. 
The electromagnetic current contains both "charge" and "polarization" contributions, and we obtain the "wormhole 
effect"— along flux lines with gapped surfaces as a consequence of the coupling. The existence of the Dirac fermion at 
surfaces when time-reversal-symmetry is present, which is the last result in Section V, becomes more understandable 
on noting that the combination of scalar and vector bosons that the bulk BF theory implies at the surface are precisely 
those required to represent a Dirac fermion via the explicit "tomographic" mapping in 2-l-lD^^!?'*. 

Some generalizations and future directions are reviewed in Section VI; one is that the BF theory description allows 
for fractional statistics not just in 2D but in 3d 15^16,25^26 _ This does not contradict the classic argument that pointlike 
objects do not have fractional statistics in 3D22 because the fractional statistics are between one pointlike and one 
stringlike object. As yet there is no microscopic realization of a "fractional" topological insulator (we discuss some 
issues in closing), but the BF approach can at least formally describe such states. Rather than try to give an overview 
of the other features of BF theory at the outset, it seems desirable to proceed with deriving it from the tt fiux response 
in the simplest case, then explain how other features of topological insulators arise in the BF theory. 



II. VORTICES WITH vr FLUX IN A 2D TOPOLOGICAL INSULATOR 



In topological insulators, an electromagnetic tt fiux has special roles both in 2D and 3D. When tt flux is threaded, 
there are always some excitations which are stable to any local T- symmetric perturbation. In 2D TIs, tt fiux 
can induce spin-charge separation with semionic statistics. In 3D TIs, tt flux creates helical ID metals. Since this 
property is topological, tt flux response can be used to distinguish TIs from trivial phases without referring to edge 
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states. However, tt flux is a strong perturbation from tfie point of view of an electron because the electron wavefunction 
gains a (—1) phase factor when it is rotated around the flux. If our goal is to depict the TI as a condensed phase, as 
in the early work on the Chern-Simons description of the quantum Hall effect, then for the condensate to consist of 
charge e objects, the phase factor must be somehow compensated if the tt flux is to be a finite-energy excitation. 

So there must be another global effect that compensates the (—1) phase factor for electrons when tt flux is threaded. 
We will show that this is the origin of a BF-type action. 

1. Review of the 2tt flux vortex 

Let's consider a vortex threaded into a 2D condensate of charge e — 1 objects (electrons). Consider the condensate 
wave function ip and a J7(l) electromagnetic gauge field A in a vortex configuration. The Lagrangian density for the 
single vortex follows = c = 1): 

L = i I (a,, - P -nV'V), V = ^{^j^^j-vy (4) 

The minimum energy configuration for ip is easy to obtain : = vexp{i9) where 9 is the angle with respect to 
X-axis in 2D. This distribution of ip has two key features. First, it minimizes the vortex potential V. Second, ip is 
single- valued because of the angle dependence exp{i0). And if we substitute this distribution of ip into equation 
we obtain just 

L = ^{d,e-A,r (5) 

The equation ([SJ fixes the field configuration of A to make the system have finite energy: Ai should be reduced into 
diO as r — >■ oo. Otherwise, the energy of the vortex is simply infinite which is not physically acceptable. Due to this 
constraint, the fiux trapped in the vortex configuration is quantized by 2n: 

J (fxeijd^Aj ^ <j) d'^ = A9 = 2tt (6) 

In general, the single-valuedness requires ip cx exjp(in9) which corresponds to an n-vortex configuration. Then this 
gives the total flux 2mT threaded in the 2D condensate. Hence, we directly see that tt flux is not allowed for the vortex 
configuration with finite energy. Thus, we need another component in the theory to deal with tt fiux vortex. 

2. TV flux vortex and BF theory 

As explained in the previous section, a finite-energy vortex configuration with tt flux is seemingly excluded for a 
charge-e condensate. This result can be traced back to the phase factor (—1) for electrons rotating around tt flux. 
Thus, if there is a way to create another phase factor (—1) for electrons, then (— 1) = 1 would allow tt flux in the 
vortex conflguration with finite energy. This is done by introducing a J7(l) statistical gauge field that couples to 
electrons. This is depicted in the figure ((T)). Now the Lagrangian ^ is replaced by the following equation. 

L = 1 I (9^ - + aJ)V P -nV-^^) V^]^(^H-v^Y (7) 

Both gauge fields A and a carry fiux tt for the simplest case. However, a can carry {2n — 1)t: flux (or equivalently, tt 
flux of a induces phase factor exp(i(2n— l)7r) = (—1)) because the single-valuedness of the wave function only requires 
the total phase change to be 2n7r in rotating around vortex. For the purpose of deriving the BF-type action, n = 1 
is enough, and we will stick to n = 1 in this section, but the generalization to n 7^ 1 is trivial, and the importance 
for n 7^ 1 will be pointed out later when we discuss the statistics ofSFtheory and the Z2-ness of TIs. With the same 
ansatz ijj = v exp{i0) as before, we get the Lagrangian density for the gauge flelds a and A: 

L^^^{d^e-A^-a^f (8) 

We can use the standard dual representation for the vortex by introducing an auxiliary U{1) fleld Then the 
Lagrangian density ([8]) can be written as the following: 



(9) 
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Now the key step of the dual representation follows: writing = 9 smooth + Ovor- The integration (or gauging away) 
of Osmooth gives the constraint 9^^'^ = 0. This can be solved by introducing another U{1) gauge field b such that 
= e'^'^^d,yb\, and this transformation is enough to show there should be a BF-type term in the Lagrangian to deal 
with TT flux vortex. After using the relation between ^ and 5, we obtain 

^ = -^/'. +£^"^5-^A(9^e.or -A^- a^), (10) 

where / is the field strength of b. Finally, we identify the vortex current as — ^s^'^^di,d\9vor, and write the 
Lagrangian in a slightly different form: 

^ = ^i^-^'- + - ^'''''(^M + <^^^)^.bx (11) 

Now we note that the first two terms in the equation (jlip are induced by the vortex configuration, and argue that 
the effective theory for the background TI should be the third term including two gauge fields a and b with non- 
dynamical external gauge field A. Another way to motivate discarding the first term is that, when we deal with 
long-range and low-energy physics, we can ignore this Maxwell-like term for the same reason that pure Chern-Simons 
theory is the effective theory for quantum hall systems: the third term in equation [TT] includes only one derivative, 
whereas the Maxwell-like term carries two derivatives and is relatively small at long length scales. With slightly 
different normalization, we obtain the effective theory for the 2D topological insulator: 

Lbf = -e'^-'^A^d^bx + a^dM) (12) 

TT 

The first term in equation (|12p is current-gauge coupling where the current is mediated by the gauge field b. The 
second term is called the BF term (because it includes the gauge field b and the field strength of the a field), and it 
contains the coupling between a and b. In fact, the BF term is topological, i.e it does not carry dynamics; rather it 
carries the information of statistics between sources of a and b. Note that the first term is actually a BF term as well, 
coupling the internal field b and the external field A. This type of mixed BF term has previously been obtained for 
the 2D topological insulator by Goryo et ali^"— . The meaning of the coefficients and the BF term will be clarified 
in the next section, and in this section we focus on consequences of (1121) . There is an interesting property for the 
equation ((T^): if we integrate out a and b, then the resulting action is 0, i.e., the bulk effective action for A is 0. So if 
we want to have an electromagnetic response due to the external gauge field A, we need another term coupled to a, 
as in the T-breaking Chern-Simons term for the quantum Hall effect. Another difference is that the BF theory with 
A above preserves P- and T- symmetry which will be clear in the next section. Looking ahead, a key feature of this 
BF theory is that it has a direct 3D generalization; for 3D, we will similarly derive the existence of a two-form gauge 
field b coupled to a and A, 

Lbf = ^e^^'^P^A^^dMp + a^^dMp) (13) 

Then the equation (fT3| inherits the properties of 2D BF theory p2|l : it is topological, the effective action of A is 
zero, and it is P- and T- symmetric. 

III. GENERAL PROPERTIES OF BF THEORY 

Our goal in this paper is to argue that BF theory is the effective theory for TIs in a way similar to CS theory for 
the quantum Hall effect. Before continuing with that argument, let us review a few more properties of BF theory. 
BF theory is a generalization of CS theory, and both are topologicali^ii^i^; in fact, the relationship between BF and 
CS theory is fairly simple in (2+l)D. This relationship is essentially the same as the microscopic relationship between 
(two copies of) the integer quantum Hall effect and the 2D TI. But there are important differences between BF and 
CS theory, such as symmetry; some similarities and differences between CS and BF theory in two and three spatial 
dimensions will be reviewed in the following. Versions of BF theory have been applied in condensed matter physics 
previously to study superconductors and topological spin liquids"'^^»^i^, but not topological insulators as far as we 
are aware. Henceforth we follow the condensed matter convention and refer to systems by their spatial dimension, 
e.g., 2D means (2-|-l)D. 
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A. 2D BF theory and CS theory 

The most general form of BF theory m 2D is composed of two U{1) gauge fields a, b and corresponding sources j, 
ji5^ The full BF action is 

Lbf = ^ef^^^a^dM - .fa^ - J^b^ (14) 

The first term in equation (jl4[) is topological and describes braiding statistics with parameter k : when the quasiparticle 
type j is rotated around the quasiparticle type J, the total wave function obtains ^ . This term has the same role as 
the term Lcs — ■^^^'^^o.^iduO-x in CS theory, and the equations of motion are also simply obtained: 

^e'^^^dM - :^£'^"^a,aA = (15) 

ZTT ZTT 

The symmetries of the gauge fields a and b are different: under T, (00,01,02) — > (00,-01,-02) and (&o,&i,&2) — >■ 
(-601^1:^2)- So the BF term in the equation (I13p is T - invariant. Armed with these, we reconsider the effective 
theory 

Lbf = -e^'''^{A^dM + a^d.b^) (16) 

Two comments about this theory are in order. First, because {Aq, Ai, A2) — > (Aq, — Ai, —A2) under T- transformation 
(the same symmetries as the a gauge field, which is why it is natural to describe the tt vortex by combining A and 
a), we see that the theory is T- invariant, i.e., the minimal coupling between A and b is T- symmetric. A acts as a 
source of b field. Secondly, our theory for topological insulators is described by A: = 2, and this makes sense because 
TT flux of a field should give the phase factor (—1) to the total wave function as we discussed tt flux vortex theory(In 
fact, k = 2 corresponds to the semionic statistics for a and b). However, this is not the only coefficient giving phase 
factor (—1) between a and b. In general, k = 2n-i should be equally good because exp(i(2n — l)7r) = (—1). So the 
most general theory with tt flux having the proper statistical effect has the form : 

Lbf = \, g^'^^a^a.bA + -e^^^A^^.^A (17) 
[2n — l)7r TT 

For any value of n, the effective action for A is still 0, so there is no (bulk) response to the external field A. Before 
getting into the details of the system's response, it should be noted that the action of 2D BF theory is just doubled 
CS theory with opposite chirality; while it is known that the 2D TI can be realized as the double copy of quantum hall 
layers with opposite chirality, the above theory is saying that this is true for some properties, such as the existence 
of edge modes, even when there is no conserved component of spin. The equation (|12p can be equally well derived 
if we start from doubled CS theory and pick appropriate combinations of the two CS fields ("charge" and "spin"). 
However, it is well known that TIs can be described only by doubled CS theories describing odd-integer quantum 
Hall effects, and we wish to emphasize here that this is really coming from the tt flux response, which is a defining 
property of TIs, and not from the assumption of decoupled spin-up and spin-down layers. 



1. 3D BF theory 



The most general form of 3D BF theory includes a one-form gauge field o^, an antisymmetric two- form gauge field 
6^^, and corresponding sources j'^, E^' ^^^i-'^^'^^i^^ . Here, the source E is an antisymmetric rank 2 tensor, so physically 
the source E represents the density of line-like objects that can be pictured as field lines for electric and magnetic 
parts of the tensor. We will later connect the magnetic and electric components of E through the equations of motion. 
The corresponding BF theory should contain information about braiding of vortices and quasiparticles. Interestingly, 
3D BF theory can be used to describe fractional statistics in 3D systems. The standard argument that pointlike 
particles cannot have fractional statistics in 3D as they do in 2D^^ is that all ways of looping one particle around 
another to detect their mutual statistics are topologically equivalent in 3D but not in 2D; as a result, 3D statistics are 
described by the permutation group, and 2D statistics are described by the much richer braid group-one consequence 
is that in 2D, there are "anyonic" particles that are neither fermionic nor bosonic. In 3D, vortices and quasiparticles 
can braid noncontractibly like particles in 2D, allowing for statistics other than fermionic and bosonic. 

Hence it is possible to impose nontrivial statistics between vortices and quasiparticles, and we will argue that this 
is what happens in 3D TIs. The starting BF theory is the following: 

Lbf = -^e^'''^Pa^^Mp - fa^ - l^'^'b^, (18) 
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As in the 2D case, the first term carries the statistics: when quasiparticle j is rotated around the vortex density E, 
the total wave function obtains a phase factor ^ . The equation of motion can be obtained by varying with respect 
to a and b : 

^e^^'^^cd^bxp = f ^e^^'^Pd^ap = S^'' (19) 

And as in 2D, k = 2 (or more generahy, k — 2n-i ) ^^^^ describe a TI for which threading with tt flux yields the 
correct phase factor (—1). The external field A acts as a source for the two- form field h if we look into the equation 
[T4l Finally, this BF theory is T- symmetric because (ao, ai, 02, 03) — (ao, — ai, —02, —03), (6oi, (^ &oi, ^ij) and 

(Aq, Ai, A2, A3) — >■ (Ao, — Ai, — A2, —A3). Thus the minimal coupling between A and & is still T- symmetric, and the 
candidate theory for 3D TIs with integer n is following : 

Lbf = V, e^-'^Papd.bxp + ^e^'^^f A^d.bxp (20) 

I\ln — iJTT ZTT 

As noted in 2D BF theory, we need a term coupled to a to have a total bulk response of external gauge field A, which 
will be obtained later. 



IV. 2D TOPOLOGICAL INSULATORS 



The physics of 2D TIs is relatively easy to understand compared to that of 3D TIs because a 2D TI can be realized 
as a bilayer integer quantum hall system with opposite T- symmetry, although this picture implicitly assumes aU{l) 
spin rotation symmetry and fails to capture the Z2 nature of the topological invariant. The role of the magnetic field 
in ordinary quantum hall systems is replaced by the spin-orbit interaction, whose symmetries require opposite integer 
quantum Hall states of up and down spins along the U{1) axis. When there is no U{1) remnant of spin-rotation 
symmetry, we are forced to interpret B as so-called time-reversal gauge field, and J as the time-reversal-pair current. 
In fact, the above Lagrangian for 2D TIs is very similar to that of quantum hall systems, and many aspects of physics 
in 2D TIs resemble those from quantum Hall states. We quickly review the topological field theory description of 
quantum hall systems. 



A. Quantum Hall systems and CS theory 

Let us consider a. v — ^ fractional quantum hall system, i.e., the lowest Landau level is 1/3 full. The system 
has elementary excitations with fractional charges : fractionalized quasiholes(= |) and quasielectrons(= — f )• These 
excitations can be made in the system if we thread magnetic flux, and the braiding statistics of those quasiparticles 
are fractional. Remarkably, these seemingly complicated properties can be simply formulated in CS theory : 

Lcs = -^e^'^^apd^a^ + fa^, f - ^^£'"'^9.^^ (21) 

47r ZTT 

Here k parametrizes the braiding of quasiparticles and flUing factor v — ^ is described k — Z. The requirement of 
gauge invariance for the total action on a manifold with boundary, with Lcs as bulk action, leads to an edge chiral 
boson field which is indeed present in quantum hall systems^i^. Moreover, if we integrate out the internal gauge 
field a, we obtain the correct effective electromagnetic Lagrangian for a state with Hall effect a^y — e^/3/i: 

LQHE = -^e^'''^A^d,Ax (22) 

Thus we can argue that CS theory is a valid effective theory for quantum hall systems due to its correct description 
of statistics of quasiparticles, correct effective theory for A, and edge theory. In the next and subsequent sections, we 
will argue BF theory also captures statistics of quasiparticles, correct effective theory for A and B, and edge theory 
of 2D TIs, i.e., BF theory should be the correct effective theory for TIs. 



B. 2D topological insulators and BF theory 



Topological insulators are gapped topologically ordered systems like quantum Hall states, and thus it is natural 
that they could be described by a topological field theory. However, BF theory, unlike the simplest form of CS theory. 



7 



involves two entities in the theory a and b with different symmetries. If we want to follow the philosophy of CS theory, 
then a and b should couple to currents of elementary excitations in the system. So it would be good to think of 
possible excitations in TIs before moving into detailed mathematics. 

To construct an analogy between CS theory and BF theory, we will stick to 2D TIs in this section. In a quantum 
hall system, we can collect the elementary excitation by threading the magnetic flux. Similarly, we can thread a flux 
of the external electromagnetic gauge field into a 2D TI, and it is known that we can have two types of excitations in 
the presence of a tt flu x^^i^" . In 2D TIs, the excitations are a charge-neutral Kramers doublet (charge and, if one 
component of spin is a good quantum number, spin= ±i) and holon(charge= ±e and spin zero). So a and b should 
couple to the sources of external gauge fields which are currents of "spinons" and "holons" ; of course the external gauge 
field coupling to the Kramers-doublet is not physical, and the other gauge field is ordinary U{1) electromagnetism. 

Note that even if no U{1) component of spin rotation symmetry like Sz is conserved microscopically, there are still 
two U{1) currents macroscopically, not just the U{1) charge current. This is reflected in the existence of two dissi- 
pationless edge modes that are asymptotically decoupled at low energy, and the generic existence of a dissipationless 
"spin Hall effect" carried by edge states in an applied electric field even if the direction and magnitude of this spin 
Hall effect are not quantized. In the BF description, there are two U{1) currents of which one couples directly to the 
external electromagnetic field; the quantized particle source of the "spin" gauge field can be visualized as one of the 
charge-neutral excitations created by insertion of a tt flux. 

We recall the previous BF theory written for 2D TIs : 

Lbf = -e^-^A^^.^A + a^dM) (23) 

TT 

From the above equation, we see that b couples to A and thus represents the field for holons. Naturally, a should 
represent spinons. More precisely, the current for a couples to a fictitious time-reversal gauge field B; the term for the 
coupling between spinons and external time-reversal gauge fields was absent in the the previous discussions because we 
only dealt with the physical EM field. To make the duality between a and b in 2D explicit, we now add -e^'^^B^d,,a\ 
to our Lagrangian density. Then the full Lagrangian for our system is obtained : 

Lbf = -e'^'^a^dM + -e^^^A^^.^A + -e^^^B^d^a^ (24) 

TT TT TT 

The equations of motion for a and b can be obtained if we vary b and a. 

e^'^dM = -e^'^d^Bx, e^'^^d^ax = -e^-'^d.Ax (25) 

From the above equations, some important features can be directly read off. First, unit spinon and holon are localized 
at the TT flux of A and B from the equation of motion. Second, this spinon and holon obey semionic statistics. 
Thus, the above BF theory ([M)) correctly describes the elementary excitations in 2D TIs. Now, consider the effective 
Lagrangian for the external fields A and B obtained by integrating out a and b. By completing squares for a and b 
and integrating out them, we obtain the following effective Lagrangian: 

Leff^-e^-'^A^d^Bx (26) 

TT 

Note that there is no bulk electromagnetic action for A alone, consistent with the absence of a charge Hall effect; 
the effective action for A and B in bulk describes a "quantum spin Hall effect" (a spin current flows in response to 
a charge fleld and vice versa) if there is a conserved current of spin in some direction. Since we can trivially rewrite 
the bulk BF theory as a doubling of CS theory in 2D, it contains two chiral bosons at the edge, consistent with 2D 
TIs. We will review the derivation of this edge theory later in Section V as preparation for the surface theory of the 
3D TI. 

Before finishing this section, it is worth mentioning briefly how the BF theory can capture the physics of the "Z2 
odd" property of 2D TIs. Consider an A^-flavored BF theory to describe N copies of 2D TI 

i=N 

L=y] i-e^-^'al^d^bl + -s^^-^'A^d^bl + -e^-^a'^d,,Bx) (27) 

1=1 

This edge of this Lagrangian can be studied from the iiT-matrix of the system. The _fC-matrix for the single-flavor 
BF theory is = diag(l, —1) when the BF term is diagonalized. This if-matrix captures that there are two one- 
dimensional edge modes with opposite chiralities. When there are N flavors, the _fC-matrix is trivially generalized as K 
= diag(l, —1, 1, —1, 1, — 1 • • • ) where the block diag(l, —1) is repeated times. This edge is unstable to T-symmetric 
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FIG. 1: The black arrows represent one-half magnetic flux quantum (a tt flux) and the black spheres represent electrons. 
Consider a system with Sz conservation for simplicity: an electron can be thought as a composite of a holon (red spheres) and 
an Sz = ±1 spinon (blue spheres), a) Whenever an electron is circulated around the half flux quantum (tt flux), the electron 
acquires a phase factor (—1) which shifts the energy levels of electrons. b)This phase factor can be cured by introducing a 
statistical gauge fleld a which is tied to the EM flux tubes; the BF term does this job. Hence the electrons see an effective 2n 
flux and behave as though there is no flux tube present, c) Spinous and holons have relative semionic statistics and thus every 
holon is affected in the same way as an electron whenever there is a spinon present. This effect generates the phenomenon that 
EM TT flux captures one spinon and tt flux of the flctitious spin-gauge fleld captures one holon, as illustrated in d) and e). 



perturbations if is even. Thus, we are led to conclude N — {2n — 1). This (2n — l)-flavor 2D TI can be fit into a 
single-component BF theory if we consider only the response to the external gauge fields. (That is, if we only look at 
the overall response to the external gauge fields, we do not actually need to retain the index i which labels different 
flavors, but can consolidate all flavors into one.) Then we are led to the simplified Lagrangian 

L = (2n - 1) X {-s'^'^^a^.d.bx + -e^^^A^d.bx + -e'^'^a^d.Bx) (28) 

TT TT TT 

This theory has a {2n— l)-times amplified spin Hall effect. This theory still has a well-defined response to the tt fluxes 
as TT flux is still the smallest T-symmetric flux. Hence, we still want to treat tt fluxes as the unit sources of a and b. 
This means we need to rescale the gauge fields a and b: {2n — l)a — > a and {2n — 1)6 b. This scaling applied to the 
Lagrangian ([28)1 gives 

L - —^-—e'^-'^^dM + -e^'^^A^dM + -e^^''''a^^,Bx (29) 

(2?! — IjTT TT TT 

which has been seen before when we considered 7r-flux responses. 

In conclusion, BF theory can describe key aspects of 2D TIs: the importance of T- symmetry, the elementary 
excitations and their statistics, and the edge theory. 



V. 3D TOPOLOGICAL INSULATORS 



3D topological insulators are described by four Z2 topological invariants, one "strong" and three "weak" . Topological 
insulators with only nontrivial weak invariants can be adiabatically connected to layering of 2D topological insulators 
and are unstable to disorder. On the other hand, the strong topological insulator is an intrinsically 3D system and 
does not have a two-dimensional limit. Moreover, the essence of the 2D topological insulator can be understood 
relatively easily by considering "quantum spin Hall" systems where a U{\) subgroup of spin rotation symmetry is 
preserved. However, if we try to preserve a U{1) spin rotation in 3D topological insulators, we end up with only 
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weak topological insulators rather than the strong topological insulator. This makes it difficult to understand the 
strong topological insulator because we need information about spins in the system but we cannot identify states by 
a remaining spin quantum number. 

BF theory supplies a route around this difficulty. Roughly speaking, the b and a fields in BF theory carry in- 
formation about spin and charge and couple to the external electromagnetic gauge field A. In this section, we will 
figure out what a and b should be on physical grounds and justify that our version of BF theory is a correct effective 
description of 3D strong topological insulators: it captures the quantized magnetoelectric effec t^"'^^ and the braiding 
of excitations in the system when the surfaces are gapped, and as mentioned in the Introduction it also allows the 
possibility of gapless fermionic surfaces. 



A. 3D topological insulators with gapped surfaces and BF theory 

As we did above for the 2D TI and BF theory, we need to identify the coupling for the internal gauge field a. The a 
field is still a one-form even in (3-|-l)D BF theory, and the external gauge field A is also a one-form. So the simplest 
coupling would be ~ daAdA, but this breaks T- symmetry of the system as T takes daAdA — —{daAdA). However, 
note that the term is a total derivative and hence can be sourced from a surface where its coefficient changes, i.e 
da A dA — d{a A dA). Thus, we are led to conclude that our theory describes TI with broken T- invariance only on 
the surface. As before, we also have T- invariant couplings between b and a, A. The resulting theory contains three 
terms. 

Lbf = ^e'^'^^Pa^dMp + ^e^^^PA^d.b^^p + Ce^-'^Pd^a^dxAp (30) 

We will find Z2-ness, the existence of only two phases as a result of T invariance, by studying the third term in 
the theory. To proceed further, we need an extra piece of information which is not encoded in the field theory itself: 
the gauge charge lattice. This is because the response to EM field is sensitive to the quantization of charges in the 
system. For convenience, we deal with a compact (3-|-l)D manifold such as T^. On such manifolds, we have no 
boundary and hence the third term might be thought to be trivial. However, this is not the case as we will consider 
T- invariance of the theory. We separate the terms in the theory (|30l) as follows : Lbuik = ^£adb + ^eAdb and 
Lsurf = CedadA where anti-symmetrization is done implicitly. The full theory is exp(z /(ibuik + ^surf))- Under the 
time-reversal symmetry operation, we obtain exp(i /(£buik ~ -^surf))- Hence, for the theory to be well-defined and 
have T- symmetry, we need to have exp(i J Lgurf) = ±1- 

Now, we impose the gauge charge lattice for the theory. First, there is a well-defined particle excitation which is 
the source for the gauge a. We quantize j° to be integral over each of the 2D tori within T*. Then from the second 
term in -Lbuik, we are led to conclude that is also quantized to be an integral multiple of e. This seemingly trivial 
imposition of the gauge charge lattice gives a nontrivial effect on Lsurf as the fluxes threaded in each 'torus' of the 
manifold should be quantized by ^ = 27r (here e = I in our units). Then Lgurf can take only the following possible 
values. 

J Lsnr f = C2ttZ X 2ttZ X 2 = CSn^Z. (31) 

Here Z indicates an unknown integer. The first 27rZ is due to the integral of dA and the second 2tt is due to da, and the 
last factor of 2 is due to the equation of motion connecting da to dA. In result, we conclude that exp(i / Lsuri) = ±1 
as C = 0, -^modj:^, and the topologically non-trivial insulator corresponds to C = 

Now, we integrate out a and b fields in the equation (pO)) with C = so that we get the effective action for the 
electromagnetic field A. When we do this integration, we obtain the magnetoelectric effect of the strong topological 
insulator—, 

S^D = ± f d^xdt^e^-'^Pd^ A.dxAp (32) 

Now we introduce the theta-angle 6 = C x 8tt which is the usual "axion electrodynamics"— of the topological 
insulator. Because of the quantization of C, we have 9^6 + 2tt. We allow d9 on the surface of the topological 
insulator, which then influences the rest of the material through the dad A term. 

With this knowledge and the full Lagrangian, we can now study other properties of the system by the equations 
of motion resulting from equation (1301) . Before getting into the details of the equations of motion for a and &, it is 
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important to notice that electromagnetic current is conveyed by both a and h. The current can be directly read off 
from the equation ([5(1)) by j^Lbf- 

J^EM = Ji: + J':- ^^'""d^hp + ^e'^-^'d^ed^a,) (33) 

So the total electromagnetic current is carried by both gauge fields a and b. Even though the two terms of this current 
appear similar, the interpretation of a and b will show that the two are independent contribution to Jem- 

It is easy to interpret the contribution from b : 6 is a two-form and thus can be thought of as related to ordinary 
particle current in the same way as the field strength in electromagnetism. Because Jb is a free-particle current in 
the TIs, we might expect Jb to appear only on the edge as the Hall current when we solve the equation of motion for 
6, and it turns out to be so. For J^, we need to be somewhat careful. The term including Ja is originally from the 
surface Lsurt- Here we treat the term as a part of bulk theory. Thus, if we only ask about the bulk current, Ja = 
trivially. But Ja need not be zero on the surface. Actually, Ja = ri-P + nxM where P is electric polarization, M is 
magnetic polarization and n is the normal vector to the normal. To see this, we revisit Lsmf for the case 9 = tt: 

Lsurt = ^e'^'^^Pd^a^dxAp = ^-P^^^d^A, = \{P - t + M - t) (34) 

OTT Z Z 

Here P^'^ is the antisymmetric polarization tensor, are usual electric, magnetic moments. We identify 

and e^^^A'P — P^'^ where P^'^ = j^e^'^^^'dxap. The polarization has an additional constraint like usual polarization 

in neutral matter: d^P^'^ = implying that V • = and V x Af = 0. From this constraint, we can see that we 
split the electromagnetic response into two pieces that cannot mix up. One part is Ja which is particle-like response, 
and the second part is P^^ which is like string or polarization response. With these in mind, we can cast our BF 
theory into the physically straightforward form 

Lbf - ^s^'^^Pa^dMp + J^A^ + ]-(P-t + fl-t) (35) 

ZTT Z 

The second term is the usual coupling of the EM field and charge current. The third term is the coupling allowed 
by symmetry between (magnetic and electric) polarizations and the EM field. The term shows up only in the 
nontrivial topological insulator. Finally, the first term gives nontrivial braiding between particle and polarizations. 
This description will give in the next section the axion electrodynamics form of EM responses of topological insulators 
with gapped surfaces that break T-symmetry once we integrate out the a and b fields. If on the other hand we did 
not break T- symmetry on the surfaces, then we will obtain a different description as seen below. 

Before finishing this section, we present the equations of motion for the fields a and b for future use. We vary a 
and b for the Lagrangian ([5(11) . The resulting equations are as follows : 

^e^'->^PdMp^±^e^'''^Pd,edxAp, ^e^^'^^Pdxap^-^e^'^^PdxAp (36) 

ZTT OTT ZTT ZTT 

Given A, we can obtain part of a and b according to these equations of motion. Effectively, A is the source for a and b. 
The most studied example for 3D TI is an EM flux tube threading the bulk and its surface. The gauge configuration 
for a is easy to capture as it describes another vortex tied to EM vortex with the same magnitude. On the other 
hand, we need to be careful for b as b has more gauge degrees of freedom. For convenience, we think of a TI with 
z < and the magnetic flux threading the surface at the origin. We stick to the static configuration of A fields. Then 
the equation for the gauge b can be split into two parts: 

J°b - ^5{x)5{y)5{z) = ^d.,e^=%,k, Jl ^ ^ ^{^ x b°'^)\ (37) 

OTT ZTT ZTT 

where we used dgb = for the second equation. These two equations imply that the end of the flux tube is a monopole 
source for ebjk, and has a curl- free configuration. This is illustrated in Fig. HI 

B. Gapped edge theory of 3D topological insulators 

In this section, we will establish the "edge theory" of the surface of a 3D TI once time-reversal-breaking perturbations 
have been added. (We use edge theory in this section for comparison to the well-known edge theories of 2D topological 
states. That the surface in our BF theory should be T- broken is visible in Lsurf-) Generically the surface is gapped 
for T- symmetry-breaking perturbations present on the surface. The gapped surface should not change the bulk 
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a) b) 




FIG. 2: A flux tube threads a 3D TI and penetrates the surface, a) The black arrow represents the magnetic flux and the black 
sphere represents the localized EM charge due to the QHE on the surface. The equations of motion for and give the 
distribution of the gauge flelds. b) The picture for a,i is simple: the gauge field generates another flux tube at the location 
of the A flux tube. The tensor gauge fleld b^^ can be broken up into two pieces 6"* and ei-jkV^ . With the flux conflguration of 
A, c) 6°' sees a monopole source at the surface and d) eij^V^ sees no source but should be curl-free according to the equation of 
motion. The picture of 6^1/ is determined only up to the gauge freedom of the fleld, which consists of adding any divergence-free 
vector fleld for 6°' and adding any curl-free vector fleld for eij^V^ . 

properties much for weak perturbations, and we expect the edge theory to reflect some properties of the bulk theory. 
In fact, 3D BF theory generates the topological term of 2D BF theory as its edge theory, plus additional terms 
reflecting the broken time-reversal symmetry. The edge 2D BF theory inherits topological properties from 3D bulk 
BF theory. In this section, we derive the gapped edge theory of 3D TI in the same way as CS theory predicts a 
gapless edge theory of quantum Hall states, and we explore the properties of the theory including the effects of an 
infinitely thin flux tube threading a 3D TIs. 

1. Review of the edge theory of CS theory 

In this section, we briefly review the derivation of the edge theory for CS theory in 2D. For simplicity, we ignore the 
coefficients and conflne ourselves to understanding the requirements imposed by the gauge invariance of the theory. 
We follow Wen's argument leading to the edge theory of fractional quantum Hall states from gauge invariance of the 
theory on a manifold with boundary^^. Let us denote the compact bulk as M and the boundary as dM . The bulk 
action for CS theory is = /j^^ a A da which is not gauge invariant. For a general gauge transformation a ^ a + dj, 
the action changes Sq ^ Sq — Jgj^j dj A a. 

Hence we introduce a scalar field (f> living on the edge and gauge transform 0-f 7 when a — a + dj. Now gauge 
invariance forces us to write the action for our system as S — So — Jgj^j d<j) A a. However, the Hamiltonian related to 
IdM dc/) A a is zero. So we need an extra term to get the usual chiral edge theory; this term reflects physics localized 
near the edge, while the zero-energy Hamiltonian reflects a topological kinetic energy part coming from the bulk. 

To introduce the edge term, we define the covariant derivative D(j) ~ d(f) — a, which is gauge invariant under any 
gauge transformation, and we write the Lagrangian S = Sq — Jqj^j d(j> A a — C Jgj^j Dip A * {D(j)) where * is the dual 
operation. The key points in deriving the edge theory can be summarized as: 1) we can restore gauge invariance by 
introducing fields at the edge, and 2) the term D(j)A*{D(j)) is not coming from the gauge invariance of the bulk theory. 
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but rather is an extra piece required for the theory to have non-trivial dynamics at the edge. This surface dependence 
is manifest in that the coefficient for this term in the edge chiral theory depends on the material; it determines the 
propagation velocity of the edge excitations. Now we use similar arguments to derive the edge theory of 3D BF 
theory. 



2. Gapped edge theory of 3D BF theory 

In this section, we derive the edge theory^^"— of 3D BF theory by requiring the total theory with surface to have 
the same gauge invariance as the interior of the 3D sample. While doing so, another important feature of the bulk BF 
term shows up; the bulk BF theory is not gauge invariant, just as above for CS theory. The surplus gauge degree of 
freedom goes to the edge and establishes the edge theory. We will initially study the case when the surface is gapped 
by a time-reversal-breaking perturbation, then return in the last part of this chapter to the question of how massless 
fermions appear when time-reversal symmetry is present. As before, we ignore the numerical coefficients at first but 
will restore the coefhcients later as the coefficients determine the statistics and Hall effect on the edge. We start with 
the edge theory without external EM gauge field A. The answer for this system can be written simply by introducing 
a new scalar field $ whose values only are significant on the boundary and writing a' = a + d^: 



bAda'~ bAa'. (38) 

M JdM 

This combined action is explicitly gauge-invariant if we define $ to transform $ — > $ — 7 when a — )• a + dj. (Note : 
$ describes degrees of freedom only on the edge since in the bulk term da = da'). It is elementary to check the gauge 
invariance when b ^ b + dS,. (We do not allow the variation with respect to b on the edge; rather we vary a and 
independently.) When we vary the nondynamical field $ to obtain a constraint, we have db — 0, i.e., b is pure gauge 
on the edge. Thus b = d( a,t least locally for the 2D edge by Poincare's Lemma and the edge action becomes simply 
2D BF theory /^^^ a A dC. 

When we include the coupling to an external gauge field A, still the same idea applies and thus we can derive the 
edge theory. If we add b A dA and 9da A dA to the bulk action, then the term b A dA requires a term b A A on the 
edge. However, 9 da A dA is itself gauge invariant, so we do not need any additional term in the edge action and we 
can move the term 9da A dA = d{9a A dA) freely to the edge for the case d9 = Q in the bulk, or vice versa. So we 
write the action for M + dM as 

{b Ada ^-b AdA + 9da A dA) ~ {b A a + b A A) (39) 

M JdM 

We can move the term 9da A dA into the edge by Gauss's Law : 

S = [ {bAda' + bAdA)- [ {b A a' + b A A ~ 9a' A dA) (40) 

Jm JdM 



The two actions ((39)) and ([401) are physically equivalent. If we integrate out a and b in equation (|39p . then we obtain 
S ~ Jj^^j 9dA A dA. On the other hand, if we integrate out a and b in equation ([3D]), we obtain S ~ J^^^ A A dA, 
identical to the case of equation ([M)) . This illustrates that we can move every physical response to the EM gauge field 
to the edge by moving the term 9da A dA. The effect of pushing this term to the edge is not only a surface Hall effect 
but also statistical, and we will discuss this effect again later. For later use we restore the important coefficients for 
the action now. We take the edge to lie in the xy plane and assume that 6 is a the pure gauge d( on the edge. Then 



S ^ f ^e^^^Pb^^dxap + ^e^^'^Pb^^dxAp + l^e^^^Pd^^a^dxAp 

- / -e^'^C^a^ap + -e^^'^PCt.duAp (41) 

JdM TI" TI" 

Some comments are worthwhile here. First, the gauge invariance of the theory is guaranteed by C/j — >■ Cm + Ca» when 
b^v b^i, + d^^i, — du^fi, and the gauge transformation of a leaves the theory invariant as a appears only as d A a. 

Secondly, when we move the term 9da A dA into the edge as 9a A dA, 9 term appearing in the edge action should 
be interpreted as 69 — 9ti — 9 vacuum, and 9 itself is not well-defined on the edge so the term proportional to da A dA 
is breaking T symmetry. Lastly, the sources of b in bulk and Q on the surface are the same. If there is a source r of 
C on the surface that enters as t^^m' then r should be a part of threading the surface (explicitly, —rii ■ 
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where n is the normal vector to the surface). Thus, there is no local excitation at the surface for 6 or C fields. With 
the full machinery, we take the edge action of the 3D bulk action as the following : 

Sedge = [ -s^^'^t^d^ap ^-E^'"' i^^d.Ap + ^s^'"' a^^d.Ap (42) 

For the bulk action of equation p2|) . we have 

S = I ^e^''^f>h^,dyap + ^e^-'^Ph^.dy^Ap, (43) 
Jm ^tt It; 

which is trivial for the A field when a and h are integrated out: in this form the entire electromagnetic response comes 
from the surface. 



3. Infinitely thin flux tube 



In this section, we calculate the charge bound to a flux tube by using the gapped edge action we derived in the 
previous section. This provides some intuition for the meaning of the fields in the effective theory. We thread a tt flux 
tube through the 3D TI adiabatically. The flux tube is inserted into a thin cylinder in 3D TIs, and the surface of the 
cylinder is itself an edge of the 3D TI. While doing so, we will use the gapped edge action because the cylinder is also 
part of the edge and thus gapped (i.e., we assume the same T-breaking perturbation continues along the cylinder). 
Unfortunately, we cannot get the gapless edge theory discussed in the next section from the gapped edge action by 
naively plugging in a tt flux configuration of A field, but we can glimpse the effects of degrees of freedom living in the 
TT flux tube of 3D TIs. 

From the action for edge (|42|). there can be at most two degrees of freedom carried by C, and a. Thus, we can guess 
there are possibly two chiral modes to the gapless cylinder surface because of time reversal properties(or no localized 
mode); indeed there are two gapless chiral modes at the infinitely thin cylinder with tt flux as expected. Now, we 
calculate the charge bound to the tt flux at the end of the flux tube. For simplicity, we imagine an infinite bulk 
with the cylindrical hole at the center and assume that the flux remains spatially separated from the surface of the 
cylinder. We take the a;-axis along the cylinder axis and the y-axis around the cylinder on the surface(We compactify 
y direction). Thus, the flux $ threaded in the cylinder is given by the integral of Ay along the y direction. 

m={>Ay{t)dy (44) 



Given this external field A, we have equations of motion for C, and a. Before deriving the equations of motion, we 
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first identify the EM current due to external EM gauge field by t^S* where S is the action for the edge: 



J'^M = fc + fa = -e^^'d.C^p + ^e^^^d^ap (45) 

^ TT OTT 

where we plugged 9 — n. The equations of motion for ^ and a are easily derived by varying with respect to a and C- 

JC = l^^^'^'^-^P^ fa = l^^'^'^-Ap (46) 
Thus, we have the following equation for the total current flowed along the cylinder. 

Q = j dtJ%M ^ J J dtdy^dtAy = ^ (47) 

Because the total flux at the end of threading process is tt, we obtain charge | localized at the one side end of the 
fiux tube and — | on the other side end of the flux. We wish to emphasize that this localized charge is exactly what is 
expected from the surface hall effect. When the surface of the TI is gapped, the surface hosts ^(half the conductance 
quantum; e = h = 1) Hall conductance, and thus we have charge | when tt fiux is threaded. One more comment is in 
order : as the surface explicitly breaks T- invariance, we can deduce that the system distinguishes tt from — tt as far 
as the fiux threads the surface. A consequence is that the charge localized at the tip of tt fiux is different from that 
of the — TT flux. However, if the flux totally lies in the bulk by making it a torus, then the system never distinguishes 
TT fiux from — TT flux. This can be easily seen from the bulk Lagrangian. 

L = ^e^'^^Pa^dMp + ^e'^^'^'A^dMp (48) 

ZiT 
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From this Lagrangian, we see that n flux is one unit of vorticity of the b gauge field. Upon encircling a unit source of 
a around the vr flux, the wave function obtains a phase factor exp(z7r) which is equivalent to exp(— ztt) where — tt flux 
is in place of vr flux. Having reproduced the expected properties for gapped surfaces, we now turn to the emergence 
of Dirac fermions at the surface of BF theory when time reversal is unbroken. 

4. Gapless edge theory of 3D BF theory 

Rather than a problem, it is a useful feature of 3D BF theory that the bulk theory by itself is not invariant under 
gauge transformations on a manifold with boundary. The gauge fields forced to live at the edge cancel out the surplus 
gauge freedom of the bulk theory. Here we will study the gapless surface of 3D BF theory and will find explicit 
expressions for a surface Dirac fermion field, constructed from the bosonic fields topologically required to exist at the 
surface of a 3D topological insulator. Combined with the ability to reproduce the expected electromagnetic coupling 
for gapped surfaces, this fermionic boundary suggests strongly that 3D BF theory is the proper topological field 
theory for the 3D strong topological insulator. 

The basic idea is that 3D BF theory contains two bosonic fields at the edgei^, one arising from and the other 
from b^^. From these two bosonic fields, it is not obvious at first sight how to construct a fermion. Fortunately, there 
are bosonized theories for (2+l)D fermions (and other higher dimensions as well'^*) which were originally constructed 
to generalize the well-known bosonization of (l-l-l)D fermionic systems. The key idea of these constructions, which can 
be made precise in terms of the "topographic" representation22i2i, lies in describing the higher-dimensional system 
as an infinite number of (1+1) dimensional 'rays'. The complicated part is in maintaining the Fermi statistics of 
particles belonging to different rays. We will find that the Dirac fermion constructed out of the two bosons predicted 
by BF theory shows the special property of the surface of 3D topological insulators that the momentum and the spin 
of electrons are locked to each other. 

In this section, we first start with the edges of 2D BF theory to clarify the structure of the surface edges of 3D 
BF theory, and move into the surface of 3D BF to obtain a Dirac fermion. While we are dealing with the surface of 
3D BF theory, we show that there are two gapless bosonic fields on the surface of 3D BF theory and proceed to the 
"tomographic" transformation of introducing rays. We can then compare the result to the bosonization of a single 
Dirac fermion. 

In order to emphasize the similarities between the 2D and 3D TIs in the BF description, we review the edge of 
the 2D BF theory first and see how the surface theory of 3D BF can be thought of as an infinite number of "rays" 
described by the edge theory of 2D BF theory. To get the edge theory of 2D BF theory, we start by solving the bulk 
equations of motion of the 2D bulk BF theory, i.e Uk = 9feA, and bk — dk^ with Coulomb gauge oq = &o = 0- Up to 
a constant, 

[ L= [ e^'^'^a^dM^ I {dtT){d,K) + {d,T){dtK) (49) 

JM JM JdM 

where M is 2D bulk, and the dynamics on its edge dM is described by the fluctuations of A and F. We see that there 
is a flrst-order edge Lagrangian from the bulk: the Hamiltonian is zero until we add additional non-universal terms 
generated by the boundary. 

The physics becomes more transparent if we take symmetric and antisymmetric combinations of A and F. Substi- 
tuting S" = (A + F)/2 and A= {k-T)/2 into the equation 

I {dtr){d^A) + {d.,r){dtA) ^ f (dtSmS) - {d^A){dtA) (50) 

JdM JdM 

The surface is assumed to generate the term 

H = l "-^{d.Af + "f{d.Sr (51) 

Where vi and V2 are the speeds of fluids and is system dependent. We identify the fluid density operators pi(x) = d^S 
and P2{x) = dxA. The equations of motions for A and S are now understood as the continuity equations for the 
fluids, 

^tp^{x,t) + {-\yv^xP^{x,t)^Q, i^l,2 (52) 

Thus, there are two chiral bosonic modes flowing in opposite directions. Moreover, we can see that the spectrum of the 
Hamiltonian (jSip is v\k\ after quantization where |fc| is the magnitude of the one-dimensional momentum, and time- 
reversal symmetry forces Vi = 1^2 • These bosonic modes can be viewed as the density excitations of a one-dimensional 
fermionic system via standard 1+1-dimensional bosonization. 
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We now wish to follow the same general prescription for the surface of 3D BF theory: we obtain bosonic fields 
and a kinetic term from the bulk action and add potential terms assumed to arise from non-universal physics at the 
surface. The result is then shown to be a bosonized version of a 2+1-dimensional free Dirac fermion satisfying the 
expected conditions for a TI surface. Starting with the BF theory with no T-breaking surface term and solving the 
bulk equations of motions for a and b with Coulomb gauge (a^ = diA, bij ~ di A Q, and ao = 6oi = 0), we obtain (up 
to an overall constant) 



L= / dtAe'^d,Cj+d,Ae'^dtCj (53) 

M JdM 

Hence the surface theory contains a scalar field A and a vector field C,, from which we need to construct fermions. The 
same surface theory for 3D BF theory has been obtained previously we are not aware that the following connection 
between BF theory and surface Dirac fermions or topological insulators has been made before. The topological part 
of the surface action for BF theories more generally is discussed in Appendix B of Ref. [s^ 
We assume that the effect of non-universal surface physics is to generate the potential terms 

H = j («i(a,A)2 ^ a2{e'''d,QY) (54) 

analogous to the potential terms at the quantum Hall edge. The constants ai and a2 are non-universal and will 
ultimately determine the velocities of excitations, which we assume to be isotropic; henceforth we set ai = q;2 = 1 so 
that the velocities of propagation are unity and we can connect to some previous field theory literature. The surface 
equation of motion from varying A is 

dts'^d^Q - V^A = (55) 

which is the total derivative of 

dts'^Q - d,A = 0. (56) 
The equation of motion from varying Q is similarly the total derivative of 

dtA~d,e'Kj^Q- (57) 

But ([55]) and ([57]) are precisely the "Bose structure" on one scalar and one vector boson that are required to construct 
massless fermionic fields via the tomographic representation; the key equation (21) of Ref. [l^ is just the above with 
-B — >■ A, — >■ and we switch to that paper's notation in the following. For self-containedness and to make this result 
less mysterious, we reproduce the key steps from that work and the explicit expression for the Dirac spinors. Note that 
the first-order kinetic energy term obtained from BF theory is not what would be obtained from, e.g., fluctuations of 
an elastic surface; a bosonization approach to the TI surface seemingly inequivalent to ours has recently been proposed 
by Vildanov^ starting from a modified hydrodynamics of an incompressible fiuid. 

We now proceed to the tomographic representation for the scalar A and vector i?., and we will quantize them to 
obtain two bosonic fields describing two incompressible fluids. The tomographic representation uses the following 
generalized delta function: 

1 1'°° 1 1'°° 

5^l^{y-n-Y)^ — d/c|fc|i/2e*''^e-*'="-"" — / dkk^'"^ cos k{y - h ■ t) (58) 

This function is used to transform the usual cartesian representations of fields into the tomographic representation. 
For example, the tomographic transformation of a scalar field 0(r) is defined as 

(j){y,n) = J drS{y — n ■ r)(j){r). (59) 

Physically, y € (— oo, oo) represents the projected distance of r from the origin along the direction n. From now on, we 
write 9 for the direction of n, i.e n = (cos(0), sin(0)). Scalar, spinor, and vector fields in the tomographic coordinates 
are denoted with an extra tilde. 

From this definition, the following two useful relations can be obtained: 

S{r -r') = ^ I dyd06{y - n ■ r)S{y - h ■ r'), (60) 
An J 
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J (frS{y ~ fi ■ r)5{y' ~ n' ■ r) ^ 5{y, y')5{n, n') + 5{y, -y')S{n, -n') (61) 

Note that the second equation manifests the connection between (n, y) and {—fi, —y) that appears below in the 
identification of the surface theory as a single Dirac fermion in condensed matter language. These are the standard 
tools for the tomographic representations of bosonization theor y^^i^'^''^^ . 

With the radiation gauge for A {W ■ A — 0), after tomographic transformation of A and B one obtains^^ 

doBiy,h) = aji(y,n) 

dyBiy,h) = d^Aiy,n), (62) 

and the longitudinal part of A is constant and can be ignored. These describe bosons propagating with nonzero 
velocity along each ray direction h. A spinor field in tomographic coordinates can be constructed for each ray as a 
normal-ordered version of 

4){y,h)^Cexp{i^/7r A{y,h) + B{y,n) , (63) 

where C is a normalization constant. The last step is to add a 2D "Klein factor" that ensures the canonical anticom- 
mutation relations of fermions are satisfied not just on one ray but between different rays; the end result is that ip is 
multiplied by an operator— 

Oft = exp (^-^ J' de' [a{n{e')) + p{h{e'))]^ (64) 

where the charges a and (3 are 



a{n) = / doA{y,n) 

J —OO 
nOO 

P{h) = / doB{y,h). (65) 



This fermion propagates with unit velocity in the +j/ direction and is the tomographic transform of a 2D massless 
Dirac fermion. A subtle point is that there is another fermion field^"^ x propagating in the —y direction, but we note 
that the Hermitian conjugate of x{—y, —h) is equivalent to ip{y, h) up to a phase factor and conclude that the above 
indeed describes just one fermion propagating in each direction; note that the surface-dependent chemical potential 
needs to be specified to determine the relative excitation energies of "electron" and "hole" propagating in a given 
direction. 

Hence the physical picture of the gapless edge of BF theory is as follows. The scalar and vector boson fields have a 
kinetic term from the bulk theory and a potential term from the surface theory, as in the Chern-Simons case. These 
combine to allow the faithful representation of Fermi fields. The details of the surface determine both the velocities on 
each ray and the filling of the Fermi states, i.e., the chemical potential. It is well understood^i^ that the remarkable 
power of bosonization for interacting systems in one dimension does not carry through straightforwardly to higher 
dimensions, essentially because the ray decomposition of the Hilbert space (a "superselection rule") is less useful 
for normal interactions. The BF theory of topological insulators predicts that interactions in the bulk still lead to 
massless Dirac fermions via the bulk-edge connection derived in this section, and additional interactions occurring 
at the surface can be more naturally treated in the emergent fermionic variables than the bosonic ones. A future 
direction that we comment on briefiy below is to consider the surfaces of fractional topological insulators in three 
dimensions. 



VI. CONCLUSIONS AND FUTURE DIRECTIONS 



We have argued that a topological field theory of BF type is the effective theory of 2D and 3D topological insulators. 
For the 2D topological insulator, this can be understood easily when one direction of spin is conserved and the system 
separates precisely into two copies of the quantum Hall effect. The BF description of the 3D topological insulator is 
our main result; perhaps the most surprising property is the emergence of surface Dirac fermions from the particular 
kinetic term of scalar and vector bosonic fields that are forced to exist by the incomplete gauge invariance of BF 
theory. 
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a) 



b) 



n 



n 





c) 



-n 




FIG. 3: a) Illustration of the tomographic representation of a distribution P = P{x,y), which is depicted as the red object. 
The tomographic projection is to project the distribution on the line along n and is described as the blue dashed density on 
the line spanned by h. b) and c) explain why {n,y) is the same as (— n, ~y), as a projection onto {n,y) is equivalent to the 
projection onto (— n, —y). 

There are several obvious generalizations and future directions that can be pursued. Considerable recent interest has 
gone into classifying topological insulators and superconductors with other symmetries than time-reversal, including 
their responses with gapped surfaces^^'''^, which are the analogues of the E-B response in the conventional topological 
insulator. Such systems may also be described by topological field theories, and it would be interesting to find the 
appropriate theories for those systems and the field-theory description of their defects and interfaces. 

We will focus in closing on one particular direction where we believe BF theory holds considerable promise. Consider 
the Chern-Simons effective theory of the integer quantum Hall effect. Multiplying the coefficient of the Chern-Simons 
term by 3 immediately gives the essential features of the fractional quantum Hall effect state previously described 
microscopically by Laughlin, including the braiding statistics and modified scaling dimension of the electron operator 
at the edge, which can be probed in tunneling experiments. 

In the same way, we can obtain an effective theory for potential fractional 3D topological insulators by modifying the 
coefficient of the BF term. It should be stated from the outset that there is not yet a microscopic parent Hamiltonian 
for a 3D fractional topological insulator. Recent constructions of effective descriptions of 3D fractional topological 
insulators have been based on "parton" ideas: the electron fractionalizes into a combination of other particles, one 
of which (a neutral spinon^^ or a fractionally charged objectA^i^^) then forms a topological insulator state. The 
BF theory presented here leads to a picture of 3D fractional topological insulators that appears somewhat different. 
Changing the coefficient of the braiding term to another of the set of possible valuesSS leads to fractional statistics 
between point-like and line-like objects. For a single pair of fields (a^, b^n)^ gapped surfaces will have surface quantum 
Hall layers of Hall conductance 1/2, 1/6, 1/10, . . ., where the first is the "integer" topological insulator studied in this 
paper and all surface Hall conductances are ambiguous by an integer. 

The most interesting aspect of such a fractional topological insulator might be its gapless surface theory. Each ID 
ray in the tomographic representation is now a pair of chiral Luttinger liquids with renormalized scaling dimension 
of the electron operator, i.e., essentially the effective theory discussed previously''^ for the edge of a 2D fractional 
topological insulator. The 2D surface state formed from the collection of such rays would be an unusual non-Fermi- 
liquid state worthy of future study. Leaving this fractional case aside, we can conclude from the results in this paper 
that BF theory can describe the universal physics of the currently extant topological insulators just as Chern-Simons 
theory captures the universal physics of quantum Hall states. 
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